Exact Matrix Product States for Quantum Hall Wave Functions 
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We show that the model wave functions used to describe the fractional quantum Hall effect have 
exact representations as matrix product states (MPS). These MPS can be implemented numerically 
in the orbital basis of both finite and infinite cylinders, which provides an efficient way of calculating 
arbitrary observables. We extend this approach to the charged excitations and numerically compute 
their Berry phases. Finally, we present an algorithm for numerically computing the real-space 
entanglement spectrum starting from an arbitrary orbital basis MPS, which allows us to study 
the scaling properties of the real-space entanglement spectra on infinite cylinders. The real-space 
entanglement spectrum obeys a scaling form dictated by the edge conformal field theory, allowing 
us to accurately extract the two entanglement velocities of the Moore-Read state. In contrast, the 
orbital space spectrum is observed to scale according to a complex set of power laws that rule out 
a similar collapse. 



I. INTRODUCTION 

The fractional quantum Hall (FQH) effects are exotic 
phases of matter that appear when interacting 2D sys- 
tems are subject to large magnetic fields. They are the 
foremost example of topologicaUy ordered phases, which 
are characterized by their pattern of long range entan- 
glement rather than by local order parameters^ In addi- 
tion to signatures such as gapless edge excitations, frac- 
tional statistics, and topological ground state degener- 
acy, recently ideas originating from quantum informa- 
tion, such as the entanglement spectrum, have become 
important tools for describing the topological order of 
these states i 2 - - — 

Given a bipartition of the system into two sub-Hilbcrt 
spaces, i5 = S)a <8> $)b, we can decompose any wave func- 
tion IvP) in terms of wave functions which live solely in A 
or B: 
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with the restriction that the 'entanglement spectrum' E a 
is real and that the 'Schmidt vectors' l^^'s form an 
orthonormal set (as do |\l/^)'s). It was argued^ and 
later numerically verified^ that when A, B are chosen 
to be regions in space, the entanglement spectrum of the 
FQH states are closely related to the conformal field the- 
ory (CFT) describing their gapless edge excitations^ In 
fact, for the 'model' wave functions such as the Laughlin-^ 
and Moore- Read (MR) states^ the states of the Schmidt 
spectrum are in exact correspondence with the states of 
the associated CFT. 

A second realm in which entanglement has come to 
play an important role is for a set of variational wave 
functions called 'matrix products states' (MPS)£ in one- 
dimension (ID) or 'tensor networks'— in higher dimen- 
sions. These are the variational states of the highly suc- 
cessful density matrix renormalization group (DMRG) 
method ) 11 ' 12 which succeeds because MPS efficiently cap- 
ture the structure of entanglement in many body wave 



functions^ The precise relationship between topological 
order and tensor network representations is a subject of 
ongoing work, but in ID at least a complete classifica- 
tion of symmetry protected topological (SPT) order for 
both gapped ID spin and fcrmion chains was recently ac- 
complished using the MPS representation of the ground 
state Jir— Given a set of sites labeled by i, each with 
local basis \rrii), an MPS \ip) is defined by a set of 'B- 
matrices', 
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The indices < a, < x to be traced over are called 
'auxiliary' indices, which we consider to be states in an 
'auxiliary Hilbert space' defined on the bonds between 
sites. With the proper normalization, the auxiliary states 
are in one to one correspondence with the entanglement 
spectrum of a cut on the bond. An important insight 
from the classification scheme is that a suitable renor- 
malization procedur e! 16 ' 18 can be defined which produces 
a representative state of the smallest possible x- F° r ex- 
ample, the x = 2 state of Affleck, Lieb, Kennedy and 
Tasaki (AKLT)^ is representative of the SPT ordered 
Haldane phased of the spin-1 Heiscnbcrg chain. 

The observed simplicity of the FQH model states' en- 
tanglement spectrum suggests they play an analogous 
role for the FQH effects as the AKLT state does for the 
Haldane phase. To pursue the analogy further, the 1 + 1D 
AKLT wave function can be written as a time ordered 
correlation function of a single '0 + ID' spin-i, which 
leads to its simple expression as an MPS whose x — 2 
auxiliary Hilbert space is a spin-i|2i The 2+1D FQH 
model wave functions can be written as the correlation 
function of a 1+1D CFT. Does it follow that the model 
FQH states have exact representations as an MPS with 
an auxiliary Hilbert space in one to one correspondence 
with the CFT, and if so, can they be implemented and 
manipulated numerically? 

In this paper we show that the model FQH wave func- 
tions and their quasiparticlc excitations indeed have ex- 
act representations as MPSs. As expected the requisite 
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structure of the model states is that their wave functions 
are the correlation functions of a 1+1D CFT, which im- 
plies essentially by definition that they are MPSs whose 
auxiliary Hilbert space is the CFT. We also explain how 
the edge excitations and ground state degeneracy arise 
in the MPS picture. 

Working on a cylinder in the Landau gauge, we can 
view the system as a ID chain of orbitals for which the 
B-matrices of Eq. ([2]) are the matrix elements of lo- 
cal operators of the CFT. We have implemented these 
MPSs numerically for both the fermionic Laughlin and 
Moore-Read states on the geometry of an infinitely long 
cylinder of circumference L, allowing us to measure ar- 
bitrary real-space correlation functions using the stan- 
dard infinite MPS algorithms. The infinite cylinder has 
a number numerical advantages, including the absence 
of boundaries, full translation invariancc and no curva- 
ture effects. Compared to the torus geometry^ - — only a 
single cut is required to study the entanglement, greatly 
simplifying the identification of the entanglement spec- 
trum. As we show in Sec. IIIII the computational com- 
plexity of the MPS representation is on the order 0(b L ) 
for b ~ 0(1). However, to achieve the same type of scal- 
ing in the traditional Hilbert space representation (say on 
a sphere£i2£~— ) would require N ~ 0(L 2 ) particles and 

T 2 

a Hilbert space dimension scaling as b . We also note 
that previously a conceptually distinct approach found 
an MPS for the Laughlin state in which there is one ma- 
trix per particle, rather than per orbital^ However, the 
construction does not generalize to other FQH states and 
again results in a complexity b N , which implies it cannot 
be implemented on the infinite cylinder geometry. 

In addition, we introduce an algorithm for calculat- 
ing the real-space entanglement spectrum of any state 
given as an MPS in the orbital basis. We first use the 
larger system sizes provided by the MPS representation 
to extract the topological entanglement entropy (TEE) 7 
using four different methods; first using the conventional 
scaling of the entanglement entropy,^ 
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,L- 1 + 0{L- 1 ), 
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and second from a similar scaling form we derive for the 
lowest entanglement energy, Eq = a E L — 7 + C(L _1 ), 
for both the orbital and real-space cuts. We find using 
Eq in the orbital cut converges most quickly, and as this 
form is applicable to other topologically ordered phases, 
it may prove useful in cases where small system sizes are a 
constraint. We are able to definitively determine 7 using 
all four methods for the q = 3, 5 and 7 Laughlin states, 
as well as the q — 2 Moore- Read state (c/. Tab. [J), which 
proved difficult in previous studies i 23 i 25 ' 26 i 28 

Finally, we perform a detailed scaling analysis of the 
spectrum for both the orbital and real-space cuts. Dur- 
ing the final preparation of this work, a recent preprint^! 
has conclusively demonstrated earlier arguments that 
the real-space entanglement spectrum of the model wave 
functions takes the form of the chiral Hamiltonian H per- 
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FIG. 1. The infinite cylinder geometry for model FQH wave 
functions. L is the circumference, x and r are the coordinates 

around and along the cylinder respectively. St = 27r . B is the 
spacing between Landau level orbitals. A real-space entan- 
glement cut between regions A, B would be made along some 
fixed r. 
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2.2568 ± 0.0003 


M < 10" 7 


Laughlin 


1.2956 ±0.0006 


2.996 


1/5 


0.672 ± 0.009 


4.96 


1/7 


0.28 ± 0.02 


6.88 


Moore-Read 1/2 


1.33 ±0.01 


0.21 ±0.01 7.77 



TABLE I. Extracted real-space entanglement velocities and 
TEE 7 for various model wave functions. v<j, is the veloc- 
ity of the chiral boson, and for the MR case, v x the veloc- 
ity of the chiral Majorana, in units of the magnetic length 
Ib = (h/qB) 1 / 2 . For the v = 1 integer quantum Hall state, 
the exact value of the velocity is known to be 4/i/tt. In the 
column for total quantum dimension T>, we present the value 
extracted via the orbital cut Eo around L = 25£b- (For v = 1 
integer case, we use the real-space cut S instead.) Refer to 
Sec. IVIIfor details on our numerical methods. 



turbed by local, irrelevant boundary operatorsi&22 This 
implies a scaling collapse of the entanglement spectrum in 
the limit L — > 00 for fixed CFT level n. The large system 
sizes available here give the first detailed demonstration 
of this principle, allowing us to extract the entanglement 
velocities for both the Majorana and U(l) modes of the 
MR state, as well as the form of the leading irrelevant 
corrections. 

In contrast, in the orbital cut each entanglement eigen- 
value scales as E a — Eq ~ L"^ which precludes the pos- 
sibility of collapsing the spectrum. This is contrary to 
earlier indications that the orbital spectrum showed the 
same linear dispersion^ though the case studied there 
was the 'conformal limit' of bosonic wave functions on 
a finite sphere. Taking the analogous conformal limit of 
the cylinder is trivial, so our results apply to the confor- 
mal limit of fermionic wave functions on an infinite cylin- 
der. Whether this is a difference between the bosonic and 
fermionic cases or due to the smaller system sizes of the 
earlier analysis is left open. 
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II. MODEL WAVEFUNCTIONS AND MATRIX 
PRODUCT STATES 

A number of gapped model wave functions, including 
those of the FQH, can be written as the correlation func- 
tions of a field theory in one lower dimension^ In the 
2+1D FQH effect, for example, the model wave func- 
tions are correlation functions of a 1+1D chiral conformal 
field theory (CFT)i^ Other examples with this structure 
include the AKLT states, the Toric code,— and certain 
BCS superconductors^ As we will illustrate in the case 
of the FQH effect, this structure implies that the state has 
an exact implementation as an MPS or a tensor network. 
The auxiliary Hilbert space of the tensor network is in 
correspondence with the Hilbert space of the associated 
lower dimensional field theory. In turn, the edge excita- 
tions and the entanglement spectrum of tensor networks 
are known to be closely related^ this relationship takes 
a particularly elegant form for the FQH effect due to the 
stringent constraints of conformal invariancc in 2Dj2I 

The simplest example is the Laughlin state on an in- 
finite cylinder, which can be written as the correlation 
function of a chiral boson (b(z)r- 



N 



^L(Zi) = ]Jsinh((z l - Zj)f 



(4) 



i<j 



exp 



N 

Vq^2 ^ Za ) ~ J d2z 



Throughout we will use z = x + it as a complex co- 
ordinate on the cylinder, where x runs around its cir- 
cumference of length L and t runs along its length, as 
illustrated in Fig. [TJ The filling fraction is v = |, the 

magnetic length is l B = , and the density of electrons 

The chiral boson cf> is a free field characterized by its 
correlation function on the plane or cylinder, 



log(z - z'), 
log sinh 



0] 



(5) 



In addition to an insertion e l ^^ Zi ^ for each electron, it is 
necessary to include a neutralizing 'background charge' 
Ob c = —iy/qp fd 2 z(j)(z). The background charge intro- 
duces some subtleties, as the branch cut in the bosonic 
propagator has a phase ambiguity equivalent to a choice 
of gauge for the electrons, which we will address at a later 
point. 

We can write a second quantized version of Eq. 
using a coherent state wave function in the vari- 
able -0 (which is a complex/ Grassmann number for 
bosons/fermions), which in the thermodynamic limit is 



# L [il>] = (0|e^ 2 ^ w * (z) |* L ) 

e jd 2 z[e^*^iP(z)-t^qp4>(z)\ 



(0) 



The notation is rather subtle as we are tying together 
two theories: the physical particles in 2+1D, with coher- 
ent state coordinate ip(z), and the path integral over the 
auxiliary space of the 1+1D chiral boson (/>, characterized 
by the correlation functions (•) j. Number conservation 
is enforced by the U(l) symmetry of the chiral boson. 

The structure of Eq. © is identical to that of a 'con- 
tinuous matrix product state' (cMPS) defined in Ref. HH, 
which we review briefly. Starting with a lattice system of 
bosons or fcrmions at sites with positions r, we first pass 
from the occupation basis {|w T )} to the coherent state 
basis by defining 



!W[Vv] = ^(Vv|iOK>. 



(7) 



Second, we note that the trace over the auxiliary states 
{a} is formally equivalent to a path integral over a ID 
system, with the B playing the role of transfer matrices. 
Anticipating the continuum limit, we assume there are 
matrices H, V in the auxiliary Hilbert space such that 
B[%j)(r)] = e^M+n^M. We can then take the con- 
tinuum limit of the MPS by analogy to the usual time- 
ordered path integral, which defines a cMPS, 



*[,p]=Tr c 



(8) 



Comparing the cMPS to the second quantized version of 
the Laughlin state, (O, we see that they are equivalent if 
we take the physical Hilbert space at each slice to be that 
of particles on a ring of circumference L, and the auxil- 
iary Hilbert space to be that of a chiral boson. In this 
case, H is precisely the Hamiltonian of the chiral boson 
(plus the background charge 0& c ), while V is a vertex op- 
erator, e 1 ^^ . This structure was also recently noted in 
Ref. [3l|, where, in the MPS language, they find the dom- 
inant eigenvector of the 'transfer matrix— of the cMPS, 
from which the entanglement Hamiltonian follows. 

While the cMPS representation is convenient from an 
analytic perspective, computationally it is desirable to 
have the discrete version expressed in the basis of lowest 
Landau level (LLL) orbitals. Letting w = e~^ z , in 
Landau gauge the orbitals can be written as 



<p n (z) oc e 



^(ir n x+|(r-T„) 2 ) 



w n e ^ { n (9) 



wher 



2-Kn p2 



l B is the guiding center for the n or 



bital. Viewing the orbitals as the sites of an infinite ID 
chain, we want to arrive at a discrete MPS as defined 
in ([2]), which requires finding the appropriate matrices 
B™ a ,. Based on the cMPS, we expect a will be in one- 
to-one correspondence with the states of the associated 
CFT. 

In order to extract the occupation number at orbital 
n, we take advantage of the fact that orbitals of the LLL 
(in the Landau gauge) are labeled by momentum. When 
acting on a many-body state in the LLL, we can then 
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replace the destruction operator ip n for orbital n with a 
contour integral around the cylinder, 



U{St) T U{8t) T U{Sr) T 



dw 
2ttI 



1 4>(w). 



(10) 



We chose to perform the integrals at r„, though with an 
appropriate change in normalization a different location 
could be chosen. 

The gauge of the cMPS, however, depends on a branch 
cut prescription for the background charge. It is conve- 
nient to choose the cut to consistently occur at some fixed 
x coordinate, such as the boundary of — L/2 < x < L/2. 
This choice of gauge does not produce the Landau gauge; 
they differ by a phase e l b . Choosing this branch cut 
prescription for <fi(z), but keeping tp n to be the destruc- 
tion operators for the Landau gauge, we find the Laughlin 
and Moore-Read states can be brought to the form 



V(w) = e V#<» (Laughlin), 

V{w) = xMe^^ (M oore-Rcad). 




CFT 



(11) 



(12a) 
(12b) 



Eq. (fTTj) looks like free time evolution governed by the 
Hamiltonian of the chiral CFT, H, punctuated by inter- 
actions at r„. As orbital ordering coincides with time 
ordering, we can pass to the Hamiltonian picture by in- 
serting resolutions of the identity 1 = ^ Q l a )( Q! l at po- 
sitions t = t„ ± e, where a label all states of the CFT. 
A resulting unit cell looks like 



(a n+1 \e- STH \a r , 



where the operator 
V 



dw 
2ni 



(13) 



(14) 



is precisely the '0 th mode' of the operator V(w), and like- 
wise 4>q is the zero-mode of the chiral boson. 

The resulting transfer operators are of two types. For 
the 'free' segments r € (r„,r„_i), the transfer operator 
is 



U(6t) 



-StH 



U(5T) a/ 3 = 8 a pe 



-8rE a 



(15) 



where a again runs over states of the CFT, with energies 
E a , and St = At the location of each site we define 

a transfer operator 



T aa > [ip, 



{a e 



00 



a') 



(16) 



(a) Structure of the orbital MPS. 

B U{%) T C/(f ) 



(b) Definition of the B-matrix. 

U(5t) T U(t 2 ) U( n ) 



\v) 



QW. 



(c)Structure of MPS with a quasihole insertion Q. 

FIG. 2. The structure of the orbital MPS. (a) U is free time 
evolution of the CFT, punctuated by perturbations T at r = 
T n . (b) The B-matrices are defined by combining U and T. 
(c) A quasiparticle is inserted by placing the matrix elements 
of the vertex operator Q in the correct time-ordered positions. 
It can be absorbed into either of the adjoining B-matrices. 



Stringing the transfer matrices together, we arrive at 
the exact MPS, 



*hM = Y[u{ST)T[<p ri 



(17) 



as illustrated in Fig. 2(a) We have suppressed the im- 



plicit summation over the CFT states a. The above is in 
'coherent state' form; to convert to the occupation basis, 
{|m)} we define the usual occupation basis B-matrices 
as 

£ S m (m!)l ip m = U{\5t) T[ip] U{\5t) (18) 



as shown in Fig. 2(b) Because the Laughlin and Moore- 



Read states are described by free CFTs, both T and U, 
and hence B, can be calculated exactly at negligible nu- 
merical cost (for the details of this calculation, we refer to 
Appendices lAl and |B|) . For an arbitrary CFT, their calcu- 
lation is more involved but nevertheless tractable using 
formulas developed for the 'truncated conformal space' 
approach to perturbed CFTs.— 

In order to obtain wave functions on a half or finite 
cylinder, one simply truncates the MPS using the vac- 
uum of the CFT as a boundary condition for the severed 
auxiliary bonds. If neutral excited states are used as 
the boundary condition, these produce the correspond- 
ing model edge excitations^ This structure is analogous 
to the spin-i degree of freedom at the boundary of an 
AKLT chain, which arises from the two choices of bound- 
ary condition for the \ = 2 MPS. 

The ground state degeneracy of the phase can also 
be understood. If the phase is d-fold degenerate, it is 
known^S that there are d-primary fields in the CFT, and 
the states of the CFT partition into d families which 'de- 
scend' from each of these primary fields. Each family is 
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invariant under the action of the electron operator V(u>), 
so it follows that a given bond can be consistently trun- 
cated to one of these d families. The d choices on the 
bond provide the d 'minimal entanglement states' 4S Ap- 
pendices |X] and [B] include a more detailed explanation. 

The 'thin torus' wave functions can also be understood 
as a limiting case of our construction! 41 ' 42 If we truncate 
the MPS by keeping only the 'highest weight' states of 
the CFT, i.e., those of the lowest energy within each 
family of the CFT, we arrive at a % = 1 MPS. We also 
note that an approximate x = 1 MPS for the Laughlin 
state was recently foundj^ in our language this can be 
understood as truncating the CFT to the states \P\ < 1. 

We now explain how the two conserved quantities 
of the LLL problem, particle number and momentum 
(sometimes called 'center of mass'), can be assigned to 
the states of the CFT. In the orbital basis we define the 
conserved quantities to be 

C = ^^qNj—p (particle number), (19a) 

j 

K = 'y~]j(c[Nj — p) (momentum), (19b) 



where j is the orbital index and we have included a filling 
factor dependent scaling (v = p/q) so that both remain 
finite in the thermodynamic limit. If a state is invariant 
under a U{1) symmetry transformation, the states of the 
Schmidt spectrum can be assigned definite charge. Con- 
sequently, the entanglement spectrum on bond h G Z + ^ 
can be labeled by pairs (C R , K n )- The states of the aux- 
iliary CFT have quantum numbers as well, in particular 
the total momentum \P\ of the CFT and the winding 
number TV of the boson (see Appendix [3] for detailed 
definitions). The pairs (N, \P\) and (Cn,Kn) are related 
by 



C n = N, 

K fl =q\P\ + \N 2 + nN, 



(20a) 
(20b) 



which explains how the previously observed offsets of the 
\P\ = levels depend on the number sector and bond 
location. 



III. CONVERGENCE PROPERTIES AND 
COMPUTATIONAL COMPLEXITY 

For numerical purposes we must truncate the MPS by 
keeping only the x most important states in the entan- 
glement spectrum. Most MPS algorithms (such as mea- 
suring correlation functions) can then be computed with 
time 0(Mx 3 ) and storage 0(x 2 ), where M is the number 
of sites involved in the measurement. We argue that to 
simulate the state at some fixed precision we must keep 
X ~ e acL / v (cL /v)~ x / 2 where c is the central charge of the 
entanglement spectrum, v is its 'entanglement velocity,' 
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FIG. 3. Numerically computing the TEE 7 for the v = 1/5 
Laughlin state (top) and the v — 1/2 Moore- Read state (bot- 
tom). 7 is extracted from both orbital (squares) and real- 
space (circles) cut, via the entanglement entropy S (filled) 
and the lowest entanglement energy Eo (empty), by perform- 
ing windowed fits to the form S(L), Eo(L) = a S E L — 7 at 
various circumferences L. The horizontal lines marks the val- 
ues of 7 = log 2? where T> 2 — 4, 5, 6 (top) and 7, 8, 9 (bottom). 
As L — > 00, the extracted value of 7 approaches their theo- 
retical values of 1 log 5 and i log 8 respectively. In the latter 
case we can see that L > 20^s is required for the TEE to be 
extracted with reasonable accuracy. (Insets) S vs. L/Ib for 
the four cases. 



and a is a non- universal constant of order 1. In contrast 
to exact diagonalization, the complexity scales exponen- 
tially only in the circumference of the cylinder, rather 
than its area. Use of the conserved quantum numbers 
drastically reduces the computational time, but does not 
alter the exponential complexity. 

Following Kitaev and Preskill's derivation of the topo- 
logical entanglement entropy,— we proceed under the as- 
sumption that the 'thermodynamic' properties of the en- 
tanglement spectrum, such as the entanglement entropy, 
take the same form as those of the auxiliary CFT. How- 
ever, there is no reason to expect the velocities that ap- 
pear will be universal, so in what follows all powers of L/v 
should be understood to have non-universal coefficients. 
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The exact status of this assumption for the orbital basis 
is somewhat unclear, because as we will show the orbital 
spectrum does not collapse to the CFT; nevertheless the 
appearance of 7 in the entropy S, the scaling form of 
the lowest eigenvalue Eo , and the collapse we find for the 
convergence of S with increased x appear to behave as 
expected. 

The density of states p(E) for a modularly invariant 
CFT is given by the 'Cardy' formula^ However, we must 
take into account the fact that only the identity sector 
of the CFT, belongs to the entanglement spectrum 
(at least for the 'minimally entangled— ground states). 
The corresponding partition function and density opera- 
tor are defined as 



p = Zr 1 Tr t e~ H ° 



(21) 
(22) 



The derivation of the Cardy formula requires a modular 
transformation, but the required partition function is not 
modularly invariant. This results in the explicit appear- 
ance of the modular S matrix, — log(S , J) = 7, where 7 is 
the topological entanglement entropy. Taking this term 
into account, the density of states is 



p(E)dE 



dE 
AE 



c)EL 



3e 



1/4 



(23) 



All the other thermodynamic properties follow from 
p(E). It is convenient to introduce the dimensionless 
variable fx, 



p 



n(c + c)EL 
3w 



1/4 



(24) 



We can calculate the partition function and entanglement 
entropy, 



:(I3) = J p(E)e-? E dE 



3v P „4 
t(c+5) LP 



dp. 



12 ^7 7T' 



7r(c + c) L 
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7 + 



(25) 
(26) 



The partition function ([23)) is evaluated via steepest de- 
scent about the saddle point p* = 



tt(c+c) L 
6v ■ 



As the 



entanglement spectrum is pi = e * /Z(l), a particular 



consequence of Eq. (|25jl is that the lowest entanglement 
level is po = e~[ aEL ~ 7 l for some non-universal Oe. A 
similar result was recently obtained in.— As illustrated 
in Fig. |31 for both the orbital and real-space cuts 7 can 
be extracted from the scaling of po with equivalent or 
better accuracy as from 5*. 
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FIG. 4. Convergence of the orbital entanglement entropy Sa 
for the q = 3 Laughlin state as the number of Virasoro levels 
kept (tia) is increased. (Inset) For various circumferences 
19£b < L < 40£b, we calculate the entanglement Sa of the 
MPS keeping only the lowest ua Virasoro levels of the CFT. 
For large enough tia, Sa converges to the exact entanglement 
entropy Soo- We expect the convergence to be controlled by 
the parameter fi — (i* oc n]/ 4 — (aL) 1 ^ 2 for some a. (Main 
figure) We plot the convergence of the entanglement entropy, 
e SA-Soo as a function of n~/ A — (aL) 1 ^ 2 , with a ~ 0.0875 giving 
a good collapse. 



The steepest descent analysis shows that the bulk of 
the probability comes from a region within 0(1) of the 
saddle point /x*. At this point the number of states x 
with E < E* is 



p(E) dE 



'2tt 



7r(c+c) L 

6 v 



(27) 



Alternately, we can define the number of CFT Virasoro 
levels required above the vacuum state, 



E*L _ (c + c) 



27TU 



24 



(28) 



To study the convergence properties, suppose we only 
keep states such that E < E\. The cutoff partition func- 
tion Z/a is 



E a 

dE p(E) e~P E 

~^(/3)ierfc(-\/2( MA -M*)), (29) 
with resulting truncation error 

£A = 1 - = \ crfc (V2(ma - P*))- (30) 

While the specific functional form may not remain uni- 
versal, it suggests that convergence is controlled by the 



dimensionless factor p\ — p,* 



,V 4 



(aL) 1 ' 2 . In the 



inset of Fig. [U we plot the convergence of the entangle- 
ment entropy 5a as a function of the number of Vira- 
soro levels n\ kept at various circumferences L. We then 
scale the data horizontally by plotting as a function of 
n 1 ^ 4 — (aL) 1 / 2 for a numerically fit value of a. Without 
any further vertical scaling, the data appears to collapse. 
This is somewhat surprising given the irregular structure 
of the orbital spectrum, but does validate the predicted 
form n a ~ L 2 . Choosing an acceptable fractional error 
for Sa, in large L limit we then conclude from the Cardy 
formula the required dimension of the MPS to simulate 
at fixed accuracy is 



X 



a ctcL/v 



(cL/v) 



-1/2 



(31) 



as claimed. 



IV. QUASI-PARTICLE EXCITATIONS 

We now discuss how to introduce quasiparticles into 
the MPS. In the 'conformal block' approach to model 
wave functions^ a quasiparticle excitation at 77 is intro- 
duced by inserting an appropriate operator Q{rf) into the 
CFT correlator, 

n^v] = (Q(v)e^i Vo{Tn) ^-^ MT ^) . (32) 

\ / CFT 

We will focus on the Laughlin and Moore-Read quasi- 
holes, for which Q is a local operator that takes a partic- 
ularly simple form, 



Q(v) = 



(Laughlin) , 



Q(rj) = a{r ] )e l4 ' (r,)l2 ^ (Moore-Read). 



(33a) 
(33b) 



Here 17(77) is the chiral part of the Ising order oper- 
ator. Quasiparticles require 'quasi-local' operators,— 
which can also be included in the MPS, but we have 
deferred their implementation. 

To incorporate the quasihole into the MPS, we first 
explicitly time order Eq. (|32[) by bringing the insertion 
2(77) between the orbitals r n _|_i > ?7 T > t„. For fermions, 
this introduces a sign for each electron in the region r > 
771- . As detailed in Appendix [Dj this sign can be written 
as s'fo/V^ where tt is conjugate to the bosonic zero-mode 
and s = ±1 for bosons and fermions respectively. 

We calculate the matrix elements of Q at r = 0, 



3/ ^Qfe)|/3), 



(34) 



and then insert Q into the 'free' evolution on the bond 
between sites n, n + 1. 



U(6r) 



where 



Tn+l 



U(T a )QU(n), r a +n = ST, (35) 
The structure of the resulting 



Vr 



MPS is illustrated in Fig. 2(c) For further details on 



calculating Q for the Laughlin and MR states we refer to 
Appendix [Dj 




FIG. 5. The real-space density p(r, x) of a q = 5 Laughlin 
state with five quasi-holes on an infinite cylinder of circum- 
ference L = 30£b- Distances are measured in units of £b- 
(The top and bottom edges are identified.) 



We have implemented the Laughlin quasiholes numer- 
ically, with a resulting density profile for a collection of 
quasiholes in the q = 5 state shown in Fig. [5] As a sim- 
ple test of the result, we can explicitly evaluate the Berry 
connection associated with the transport of one quasihole 
around another, 



6 = <pdA= cbdr)(ri\{-id v )\r)). 



(36) 



We keep one quasihole fixed at 77 = 0, while a second 
follows a discretized path r)i chosen to wind around the 
other, which defines a discretized connection e lAij = 
(Vi\Vj)- We then integrate the connection after sub- 
tracting out a similar phase in the absence of the sec- 
ond particle. Calculating the inner product between 
two matrix product states can be computed with com- 
plexity 0(A/£ 2 3 x 3 ), where A is the area of the region 
enclosing the quasiholes in question. Working on an 
L = 16£b cylinder and ensuring the quasiparticles re- 
main at least a distance of 8£b apart, we find a statisti- 
cal angle 6 q= z = 2.0992, compared to the prediction of 
^ = 2.0944. The computation takes about 1 minute. 

While the result is already well established for the 
Laughlin states^ it would be worthwhile to explic- 
itly calculate the non-abelian Berry connection for the 
Moore-Read quasiparticles. As we have computed 
the form of the Q-matrices, we believe this would be 
tractable. 



V. REAL-SPACE ENTANGLEMENT 
SPECTRUM 

Finally, we present an algorithm for computing the 
real-space entanglement spectrum (RSES) of quantum 
Hall states on both finite and infinite cylinders. In con- 
trast to the orbital cut,— which cuts the system into two 
sets of LLL orbitals, in the real-space cut we cut the 
system into two regions of physical space. Previously 
this was accomplished analytically for the free v = 1 



47,48 



numerically using Monte Carlo**, and large 



scale SVD of explicit wave functions 



2«.2!l 



Our technique 
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(a)Split operation. 
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(b)Swap operation. 

FIG. 6. The splitting and swapping procedure performed on 
an MPS. The two operations allow us to compute the real- 
space entanglement cut from an orbital MPS. 



center of the set of sites we have split), we sort the MPS 
through a series of 0(M 2 ) swapping procedures, bringing 
all indices associated to the left region to the left of the 
cut, and likewise for the right. As illustrated in Fig. |6(b)] 
for each swap we form a two-site wave function, permute 
the right and left legs to bring them to the desired order, 
and then split the wave function using SVD to obtain 
a new pair of B-matrices. Again, the canonical form 
of the MPS is preserved during this procedure, so after 
performing the required swaps the bond designated as 
the cut gives the real-space Schmidt decomposition. 

Depending on the initial bond dimension \- it m &y be 
necessary to truncate the new £?-matrices by keeping only 
the largest singular values of the SVD. It appears that 
the low lying states are not affected by truncation of the 
highest lying states, but the convergence with increased 
X should be checked on a case by case basis. 



is not specific to the model wave functions, and provides 
a means for computing the RSES of non-model states 
calculated from DMRG. As the scaling form of the en- 
tanglement spectrum only appears for the real-space cut, 
this may prove an important diagnostic for non-model 
states. For simplicity, we assume a wave function in the 
LLL, and consider an entanglement cut running around 
the cylinder at r = r c . 

The first step of the algorithm is to 'split' each orbital 
ipn(z) [Eq. ©] into components (p n R/L supported on the 
right and left of a cut at r c , 

tp n (z) = 0(t - T c )ip n (z) + 6(t c - r)ip n (z) 

= g n L¥>nL(z) + g nR tp nR (z). (37) 

We normalize the split orbitals by factors g^ /r such that 
{fnL, frill} remains an orthonormal basis. Expressing 
the state in terms of the 'split' basis amount to appending 
isometrics / onto the _B-matriccs, 



B. 



a/3 



r>m jkl jkl _ 




k I r 
9L 9r Ok+l,r 



(38) 



where the factors g implicitly depend on the orbital loca- 
tion. Because the orbitals ip n arc exponentially localized 

about r« = 2 \ B n, we can work at some fixed accuracy 
by splitting only the M ~ O^j^) orbitals nearest to the 
cut. In practice, we find M = 1.5L/£b is sufficient to 



obtain a converged spectrum. As illustrated in Fig. 6(a) 



the affected B matrices are then split using a singular 
value decomposition (SVD) equivalent to the truncation 
step of time evolving block decimation (TEBD)<££ As 
with TEBD, the splitting step preserves the 'canonical' 
form of the MPS, implying that the bipartition about the 
new bond is a Schmidt decomposition. 

After the splitting step we have added M B-matrices to 
the chain, with orbitals alternating between the left and 
right sides of the cut. Choosing some particular bond to 
represent the location of the cut (usually the bond at the 



VI. ENTANGLEMENT SPECTRUM: ORBITAL 
CUT VS. REAL-SPACE CUT 



In this section we study the scaling form of the RSES, 
then contrast it to the orbital spectrum. As illustrated 
for the q — 3 Laughlin state in Figs. [7J and for the MR 
state in Figs. QTJ the orbital and real-space cuts agree in 
their counting, which is that of the CFT, but differ in the 
scaling of the energy levels Ei present in the spectrum. 

Kitaev and Prcskill^ first noted that the known univer- 
sal features of topological entanglement entropy could be 
explained if the energies of the entanglement spectrum 
coincided with those of the chiral CFT. A physical argu- 
ment was later provided in Ref. @. Recall that the states 
of the CFT are grouped into 'families' associated with 



each primary field 



39 



in this context one per degen- 



erate ground state on a cylinder, and let P$ h denote a 
projection operator onto the corresponding family. The 
basic conclusion of Ref. HI was that the reduced real-space 
density matrix of a topological state with gaplcss chiral 
edge modes takes the form 



PL 



h 



PhP$ h e 



-vH+0(kl E 



(39) 



Here H is the Hamiltonian of the CFT, which we will 
take to have velocity 1, so an 'entanglement velocity' is 
included as a factor v. H is perturbed by more irrelevant 
boundary operators of order (Mb)" for v > 1. The coef- 
ficients ph depend on the degenerate ground state being 
considered. During the final preparation of this work, 
this scaling form was put on firm footing for the model 
FQH states^i 

While the irrelevant operators generally introduce 'in- 
teraction terms' to the entanglement Hamiltonian, to il- 
lustrate the expected behavior we consider the simplest 
type of correction, a dispersive term. For the Laughlin 
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FIG. 7. Entanglement spectra (in the neutral charge sector) 
of the orbital (left) and real-space (right) cut of the q — 3 
Laughlin state at L — 32£b- The energies E are plot against 
their momentum in units of Afc = 2 j-. Both spectra have the 
counting 1, 1, 2, 3, 5, 7, 11, etc., consistent with that of a 
chiral boson CFT. However, the energy levels E have vastly 
different quantitative behaviors in the two cases, which we 
investigate in Fig. [8] The dashed line on the right is of the 
form ve(k) = vk[l + 112k 2 + 114k 4 ], with 112,^4 fit from the 
highest level of each sector, which we associate to the state 
a\ JO). The fits appear to rule out a similar term tti, but larger 
sizes and a treatment of the 'interactions' would be required 
to rule at U3 if it is indeed absent. 



state this takes the form 



E a =E a -Eo~ v 



2J e{k n )a\ l a n + 



,n>0 



2tt 



, (40) 



e(k) = k[l + u 2 k 2 + u 4 /c 4 H ], k n = — n. (41) 

E 

which accounts for the 'branches' apparent in the real- 
space spectrum (Fig.[7]right), each of which is associated 
with the presence of a new mode (V n . The dispersion re- 
lation can be fit from the heights of these branches. Note 
that only odd powers of k can appear in the dispersion of 
a chiral boson. In general, if the irrelevant perturbations 
descend from the identity boundary operator, only odd 
powers in k should appear. 

In order to accurately extract the entanglement veloc- 
ity, we consider the scaling of the shifted spectrum 
with increased L. Based on these scaling ideas, a state 
with momentum k = ^-{n^ + n x ) should have an energy 



EgL 

2ir 



v^A^h + n0 ) + v x (A xh + n x ) 



t2a,E 



t4 a L 4 + 



(42) 



where and n x are integers corresponding to the mo- 
menta of the U(l) and Majorana sectors. The offsets 
A0/ x /j are the scaling dimensions of the highest weight 
state in the sector, which depends on the bond and num- 
ber sector in question (for the Laughlin states, it is ^-). 
For the MR case we have included a detailed exposition 
of this structure in Appendix [Cl 
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(a)Real-space cut: plot of EL/2n vs. — 1/L 2 . 
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(b)Orbital cut: log-log plot of EL/2n vs. L. 

FIG. 8. The relative entanglement energies for the real-space 
cut and orbital cut of the Laughlin state at v — 1/3. (Data 
shown for the charge neutral sector, L shown in units of £b-) 
The states at different momenta are distinguished by their 
colored symbols, (a) For the real-space cut, we plot ^\ the 
relative entanglement energy (relative to Eo) times the cylin- 
der radius, as a function of — 1/L 2 . The energies are extrapo- 
lated to L = 00 using a quadratic fit to 4- t2L~ 2 + tiL -4 , and 
land on multiples of the entanglement velocity v ~ 1.2956. 
(b) For the orbital cut, we show on a log-log plot, show- 
ing that the data has a linear behavior with negative slopes. 
The lines shown results from a linear fit to the last few data 
points. This demonstrates a nonlinear relation E a oc L~ C ' a 
with ( a > 1. 



Focusing on the identity sector = of the real-space 
q = 3 Laughlin cut, Fig. 8(a) tracks the scaled relative 
entanglement energy levels as a function of L~ 2 , ex- 
trapolating their value as the circumference approaches 
infinity. As indicated by the right-most tics of the fig- 
ure, ^ approaches n^v^ for large L, where ty, ps 1.2956. 
The data clearly confirms that the real-space entangle- 
ment spectrum approaches a linear dispersion with fixed 
velocity, and the success of the fit justifies the absence 
of L _1 and L~ 3 pcrturbativc terms. We have tabulated 
the velocities for the q = 1, 3, 5 and 7 in table HI The 



relation v q —i/v q —3 ~ \/3 noted previously^ appears not 



10 



6 - 



^1 t= 

1[£] |(M 



4 - 




2 - 



Q I i i i ■ I i i i 1 1 

20 25 30 



FIG. 9. Extrapolating the entanglement energies EL/2n for 
the Moore- Read state at v = 1/2, using a cut associated with 
counting 1, 1, 3, 5, 10, etc. for N — charge sector (c/. 
App. IC|). (Left) Here we show that the energies for the first 
three momenta extrapolate to integral combinations of and 
v x , the velocities of the chiral boson and Majorana mode re- 
spectively. The major tics on the vertical axis labels multiples 
of the minor tics label combinations n^v^, + n x v x for in- 
tegers and n x (and £b set to 1). States with momentum 
4^ extrapolate near, but not exactly, to the theoretical pre- 
diction, which we attribute to smallness of the system sizes. 
The superimposed lines are quadratic fits over the largest few 
circumferences, extrapolating to give Vtf, « 1.33, v x m 0.21. 
(Right) The theoretical placement of the energy levels for the 
state. Here the boson counting (1,1,2,3,5) and the Majorana 
counting (1,0,1,1,2) are apparent. See App. [C] for detailed 
explanation of the counting in this plot and data for other 
charge sectors. 



to continue to higher q. 

This same technique can be used for more complicated 
wave functions such as the Moore-Read state, as shown 
in Fig. [9] We extrapolate the velocities of both the charge 
and neutral modes to be ~ 1.33 and v x sa 0.21 respec- 
tively. We note that the extrapolation is only possible for 
sufficiently large circumferences L > 20£b, which is well 
within reach using the MPS representation of the wave 
function. 

Figure [8(b)| shows that in the orbital-cut, E does not 
extrapolate to the CFT linear dispersion. Rather, they 
appear to follow power law decays E a ~ L~^ a with dif- 
ferent £ a for each state a. For example, the fit for k = =j- 
gives £ « 3.0, while C ~ 2.3, 2.1 for the two set of states 
at k = 2^j£. (In the real-space case, C = 1 f° r an the 
levels.) Unfortunately, the range of data available is in- 
sufficient to draw any conclusions. 

Finally, we note that one can extract the topological 
entanglement entropy in either types of cut. This was 
shown in Fig. [3] where we used both the entropy S and 
the zero momentum state Eq as a function of L. For each 
L we perform a windowed fit; presenting the intercept of 



the best line fit through the neighboring points. While 
it is possible to extract 7 from any of the four computed 
quantities, we can see that the real-space cut is less oscil- 
latory than the orbital cut. At the same time, using the 
orbital cut Eq seems to give a much better convergence 
of 7 than any of the other methods, i.e., the system size 
L required to computed T> = e 7 via the orbital Eq to 
accuracy ±0.5 is the smallest. Table Q] lists the entangle- 
ment velocities and TEE for various Laughlin states and 
the q = 2 MR state, as well as the velocities extracted 
via the method used in Figs. 8(a) and [HI 



VII. CONCLUSION 

We have shown how the CFT structure in model FQH 
wave functions enable us to represent them as matrix 
product states. These MPSs can be evaluated numeri- 
cally on an infinite cylinder; the distinct advantages of 
this geometry, as well as the efficiency of the MPS, al- 
low us to study in detail the scaling properties of the 
Moore-Read entanglement spectrum, including a defini- 
tive identification of the U(l) and Majorana modes and 
their velocities. 

There are several future directions. The MPS rep- 
resentation is well suited for studying the screening 
properties of the states as well as their Berry connec- 
tions, so it would be valuable to numerically implement 
the MR quasiholcs in order to verify various screen- 
ing arguments^ As we have noted, our construction 
also generalizes to other topological phases whose model 
states can be expressed as a correlation function of a 
lower dimensional field theory. The resulting picture is 
strikingly similar to the 'entanglement renormalization' 
classification of ID phases exemplified in the AKLT state. 
In particular, it would appear that the fixed points of 
the entanglement renormalization scheme may be inter- 
preted as some form of fixed point for the auxiliary field 
theory when expressed as a tensor network - for topo- 
logical phases, a massless fixed point, while for trivial 
phases, a massive fixed point. Making this connection 
precise would be an intriguing development. 

We would like to acknowledge helpful conversations 
with Joel E. Moore, Tarun Grover, Frank Pollmann, 
and Sid Parameswaran, as well as support from NSF 
GRFP Grant DGE 1106400 (MZ) and NSF DMR- 
0804413 (RM). 
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Appendix A: Evaluation of B-matrices for Laughlin 

states 

Here we provide more detail on the precise form of the 
Laughlin MPS and its numerical implementation. The 
mode expansion of the chiral boson is 



4>{w) 



bo, no] 



a n + 4>o + — log(tu), (Al) 
n\ % 



+m ■ 



The field is composed of the fluctuating part <p'(w) and 
the 'zero mode', 



(to) = 4>'{w) +4>o + — \og(w) 



(A2) 



The states of the fluctuating sector can be labeled by 
occupation numbers, which we denote by a string of pos- 
itive integers P. For example, |0) denotes the ground 
state, |221) = ^a\a\a\\ti), etc - Wc define \P\ to be the 
total momentum of |P), given by the sum of the inte- 
gers. The states of the zero-mode sector are labeled by 
the eigenvalues of ttq. For convenience, we define 'charge' 
by N = y/qno, chosen such that the electron has charge 
q. The states of the zero- mode sector are labeled by \N), 
so the full CFT is then spanned by \P,N). 
Treating first the 'free' evolution U, we find 



L 



\P\ + liV 2 
2q 



U(5t)p !N ;P',n' = 5p,pi5 N , N >e 



, e -(^)>l+^ 2 ] 



(A3) 



(A4) 



Now we calculate the on-site term T, first by convert- 
ing from the coherent state form T[tJj] to the occupation 
basis, T m : 



Hence A n is simply the nth coefficient of a Taylor ex- 
pansion in w. The matrices A are non-zero only for 
P — P' = —(N + N')/2, due to momentum conserva- 
tion. Numerically, we impose a cutoff A such that we 
only keep states \P) with \P\ < A, which allows us to 
evaluate A for only a finite number of states. The time 
to compute A is proportional to its number of entries, so 
the construction of the MPS is an insignificant part of 
the computational cost (i.e., compared to matrix multi- 
plication). Combining the zero- mode and fluctuations, 

_ N + N f 

(P,N\V Q \P',N')=A pp > 5 N - N ,, q . (A10) 

Finally, the sandwiching background charge contributes 
e ~7f _ Sp.p'Sn-n' .-i to each site. 

Focusing on the case of fermions where there is at most 
one particle per orbital, 



Tp,N:P',N' — $P,P'fiN-N>,-l 
w+jv' 

Tp N\P' ,N' = App/ 2 <5/V — N'.q— 1 



(unoccupied), 

(Alia) 

(occupied) . 

(Allb) 



For case of bosons, the higher occupation states involve 
products of the A's. 

The q fold ground state degeneracy of the Laughlin 
states can be seen by noting that on a particular bond, 
e 2mN/q j g a cons tant, and can be chosen to take q values. 
In the small L limit, these q states evolve into the 'thin 
torus' wave functions^ in which we keep only the \P\ = 
states of the CFT. 



Appendix B: Evaluation of B-matrices for 
Moore-Read state 



T' 



^2T m (m\) 3/2 ip m . 

m 

i 



"ml 



(A5) 
(A6) 



We next compute the matrix elements of the vertex op- 
erator, 

(P,N\V \P',N') 

= (P, N\e i ^' l '' ( ' w ^ +i ^^ a+Nlos( - w ^ IP', N'}. (A7) 

The zero-mode part depends only on N, 

( N \ e ^4>o+Nlog( w ) \ N ,j = SN _ N , q w N+N'/2_ ( A g) 

The fluctuating part depends only on the oscillators \P), 
so we define 



An 

A P,P' 



ml—.w-^e^'^^P 1 ). (A9) 



The CFT associated with the Moore-Read state is a 
tensor product of a chiral boson </> and a Majorana mode 
X- We first give a brief review of the structure of the 
chiral Majorana CFT on a cylinder^ The states form 
four sectors according to their boundary condition (be), 
(periodic 'P' or antiperiodic 'AP') and number parity 
(even '+1' or odd '— 1'). We denote the lowest energy 
states of these four sectors by l |l)' for AP/1, l \x)' for 
AP/-1, l \a)' for P/l and '|u)' for P/-1. In the peri- 
odic sector the Majorana has modes Xn ■ n € Z, while 
in the anti-periodic sector it has modes Xm :in£2 + i 
The states of the P / AP sectors can be obtained by act- 
ing with the P/AP modes X-m on I* 7 )/!!) respectively. 
Within a given sector, the states can then be labeled by 
a string of numbers P x ; they are either integers or half- 
integers depending on the be, and do not repeat because 
of the fermionic statistics. Letting \P X \ denote the total 
momentum of the Majorana, 



IP 



2n 



(Bl) 
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where A = {0, yg} for the AP and P sectors respec- 
tively, though A can be ignored as it only changes the 
normalization of the state. 

For simplicity we consider the case in which the physi- 
cal particles are fermions, which we take to have periodic 
boundary conditions. The operator e 1 ^^^ x( z ) must be 
periodic in z at the location of the Landau orbitals r ra 
(our choice of gauge has a twist boundary condition in 
between). This introduces a constraint between the zero 
mode of the boson, N — ^/q^o, and the boundary con- 
dition of the Majorana: for P, we must have N G Z, 
while for AP, we must have N E Z + |. The boundary 
conditions will correspond to different degenerate ground 
states, with 4 states of type AP and 2 of type P, for a to- 
tal of 6 on the infinite cylinder (for a torus, the P sector 
acquires an additional 2 states depending on the parity 
of the electron number) 42 

The total energy of the combined CFT is 
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FIG. 10. The thin torus orbital wave functions for the (a) 
antiperiodic and (b) periodic sectors. Each site corresponds 
to an orbital, which is either filled (black) or empty (white). 
Each bond has only a single state of the CFT, which we de- 
compose into the Majorana part, shown above the bond, and 
zero mode of the boson N, shown below the bond. 
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2n 
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IP 



\Pa 
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(B2) 



„-<5-rP 



in the chain changes the CFT labels we assign to the 
Schmidt vectors, but does not actually change the phys- 
ical state. In the small L limit, these 3q states evolve 
into 'thin torus' wave functions 42 In this limit, we re- 
strict the Majorana CFT to the states |x), W)i I/- 4 )}; 
is diagonal if we work in the oc- and fix \P$\ = for the boson - the 8 resulting states 



1 

Yq 



with \P$\ and N arising for the boson. As for the Laugh- 
lin state, U 



cupation basis. Constructing the T matrices proceeds as 
for the Laughlin case, but we must include the Majorana 
sector in the computation of V = § ^w^ 1 x^e 1 ^^. 
Letting 



(Px\Xm\P') 



(B3) 



denote the matrix elements of the Majorana operators, 
the required matrix element is 



{P^,P x ,N\Vv\PLP'N' 



?N-N',q 



, (B4) 



with A defined as for the Laughlin case. 

For fermions, where there is at most one particle per 
orbital, 



are precisely the 'highest weight' states of the CFT, as 
illustrated in Fig. [TOl 



Appendix C: The counting of the Moore- Read state 

As explained in the last section, the chiral Majorana 
CFT may be separated into four sectors, by periodicity 
of the boundary as well as the particle number parity. 

In the periodic sectors a and /x, the excitations have 
momenta which are integral multiples of Afc = ^ , hence 
the counting of level n is the number of partitions of n 
into an even/odd number of distinct non- negative inte- 
gers. The number of states at momenta 0, Afc,2Afc,... 
are as follows^ 

fi,<r: 1,1,1,2,2,3,4,5,6,8,10,12,15,18,... (Cla) 



t (Vp x ,w),(P',p^at') 

T {P^,P X ,N),(P',P'N') 



(B5a) 
(unoccupied), 

(B5b) 



= /.Xp m p' A Pj . P ' 2 &N-N<,q-i (occupied). 

ra 

For the case of bosons, the higher occupation states in- 
volve products of the Vo's. 

The 3(7 fold ground state degeneracy of the MR states 
can be seen by first choosing a be sector, P or AP. In 
the AP sector, (— \~) F e lvN / q is a constant (where F is the 
Majorana number), which takes 2q values, each leading 
to a distinct state. It is also constant in the P sector, but 
here the 2q values only lead to q distinct states. This is 
because inserting the Majorana zero-mode xo somewhere 



(Because the presence of the zero-momentum mode, the 
counting of the two P sectors are identical.) 

In the antiperiodic sectors the excitations have mo- 
menta which are integer-plus-half multiples of Afc, or in 
other words, twice the momentum is always an odd multi- 
ple of Afc. Hence in the 1 sector the counting of level n is 
given by the partitions of 2n into positive odd integers^ 

1: 1,0,1,1,2,2,3,3,5,5,7,8,11,12,... (Clb) 

The same definition also hold for the x sector, with count- 
ing as follows,— 



X 



1,1,1,1,2,2,3,4,5,6,8,9,12,... (Clc) 



Note that since there are an odd number of excitions, 
the lowest energy state is \x) = X1/2II) with momentum 
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FIG. 11. The entanglement spectra of the q — 2 Moore- 
Read state at L = 25£b, with the orbital (left) and real-space 
(right) cut, in the N = (top) and N = 2 (bottom) charge 
sectors. In the orbital case, the cut takes place on the bond 
with 1) Majorana and N = boson state in the thin torus 
limit (c/. Fig. I10[) . The real-space cut takes place at the r 
centered on that bond. For N = 0, the counting of the states 
are 1, 1, 3, 5, 10, 16, etc, while for N = 2, the counting are 1, 
2, 4, 7, 13, etc. 
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FIG. 12. Extrapolating the entanglement energies EL/2n for 
the Moore- Read state at v = 1/2. The real-space cut is phys- 
ically centered on the bond with 1 sector (in the thin torus 
limit), but the data is shown in the N = 2 charge sector, and 
hence the counting matches that of the x sector. 

(Left) The energy states at momenta 3|j, 5|j, 7|j are la- 
beled by their shape, with their values extrapolated to L = oo. 
(Right) The theoretical energy levels of the MR state for the 
X sector. (Upper right inset) The least-momentum state in 
each charge sector N at L — 25£b, with the sector of interest 
marked red. 



The lowest states extrapolates to v^, + 



consistent with 



Eq. (|42[) (using A^h = -f- = 1). (States with momenta 7^ 
do not extrapolate to their theoretical values due to insuffi- 
cient system sizes.) Contrast this plot to the N = sector of 
Fig.U 



^Afc. The counting in this sector corresponds to the 
number of states at momenta ^Afc, |Afe, |Afc, .... 

Combined with the chiral boson, the counting of the 
Moore- Read edge spectra ar o 5758 

1: 1,1,3,5,10,16,28,43,70,..., (C2a) 
X- 1,2,4,7,13,21,35,55,86,..., (C2b) 
^,cr: 1,2,4,8,14,24,40,64,100,.... (C2c) 

(Again, in the x sector, the momenta are shifted by 
hAk.) Figure [TT] shows the orbital and real-space cut 
of the q = 2 MR state giving the 1 sector. 

Notice that in the AP case, the Majorana sector al- 
ternates between 1 and x sectors whenever an orbital is 
filled (see Fig. ITU]) . Hence the entanglement spectrum 
with different charges would also alternate between the 
countings (|C2a|) and (|C2b|) . shown clearly in Fig.[TT] Fig- 
ure [JJ] shows the relative entanglement energies of q = 2 
MR state at the N = ±q charge sectors, contrast this to 
the N = sector of Fig. M 



Appendix D: Evaluation of Q-matrices for Laughlin 
and Moore-Read quasiholes 

Evaluation of the Q-matrices for the Laughlin state 
can be done in a similar manner to the bulk i?-matrices, 
but omitting the contour integration: 

Qp.N;P>m> = Sn-N',i (sw) N+N '^(P\e^'^^\P') 

(Dl) 

where w = e~~ Vx and s = ±1 for bosons or fermions 
respectively. Note that the momentum is no longer con- 
served. 

The Moore-Read case is more complex. In the context 
of the Majorana CFT, the Ising order and disorder fields 
a, fi are 'twist' fields, interpolating between AP and P 
periodic bc's. We take the point of view that the fields a 
and \i have fixed fcrmion parity +1 and —1 respectively. 
The resulting fusion rules are 

MM = \pM = [x][x] = [a], (D2) 
MM = i, MM = M, MM = [4 (D3) 
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In this approach, there are two possible quasihole inser- 
tions, cr(?7)e 1 */^, and fi^e" 1 ^^. As xo°" ~ M; this is a 
direct realization of the picture in which each vortex has 
a Majorana zero-mode. The non-trivial vector space of 
quasihole excitations arises from the freedom of choosing 
a or [i, subject to the constraint that they fuse properly 
to the vacuum. 

Fortunately the techniques for evaluating matrix ele- 
ments of the type 

have already been developed in the 'truncated-fermionic- 
space- approach' to the perturbed Ising CFTi£2 

Consider, for example, the AP to P case. Arbitrary 
states can be built by acting with the modes X-m so 
without loss of generality we consider the matrix element 

H II xmMn) II x-nM)- (D5) 

{ mi eP x } {ni£P^} 



The chief technical result of Ref. [H Eqs. 2.9-2.13 is that 
there exists an easily computed matrix C(rj) such that 

H • ■■XmO-{rf) ■■ ■ |1) = (<t\ ■ ■ ■ o-{ri)C mn {rj)Xn • ■ ■ 

(D6) 

After commuting all \ across the insertion, the Majo- 
ranas are brought to normal ordered form, reducing the 
problem to Wick contractions and the matrix elements 

(a\a\l) = C a<rl = 1, (D7) 
= C^ x = 1, etc. (D8) 

We have not as of yet implement the MR quasiholes nu- 
merically, which would be a worthwhile check given the 
subtleties of this case. 
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